Wide binaries as a critical test of Classical Gravity 
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Modified gravity scenarios where a change of regime appears at acceleration scales a < ao have 
been proposed. Since for IMq systems the acceleration drops below ao at scales of around 7000 AU, 
a statistical survey of wide binaries with relative velocities and separations reaching 10'* AU and 
beyond should prove useful to the above debate. We apply the proposed test to the best currently 
available data. Results show a constant upper limit to the relative velocities in wide binaries which 
is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation 
curves in the same a < ao acceleration regime. Our results are suggestive of a breakdown of Kepler's 
third law beyond a ~ ao scales, in accordance with generic predictions of modified gravity theories 
designed not to require any dark matter at galactic scales and beyond. 
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I. INTRODUCTION 

Over the past few years the dominant explanation for 
the large mass to light ratios inferred for galactic and 
meta-galactic systems, that these are embedded within 
massive dark matter halos, has begun to be challenged. 
Direct detection of the dark matter particles, in spite 
of decades of extensive and dedicated searches, remains 
lacking. This has led some to interpret the velocity dis- 
persion measurements of stars in the local dSph galaxies 
(e.g. 0, [3) I the extended and flat rotation curves of 
spiral galaxies (e.g. [2^, [HI), the large velocity disper- 
sions of galaxies in clusters (e.g. [11]), stellar dynamics 
in elliptical galaxies (e.g. [3J|), the gravitational lens- 
ing due to massive galaxies (e.g. [ll], [1]), and even the 
cosmologically inferred matter content for the universe 
through CMB and structure formation physics (e.g. (soj . 
0, [3), not as indirect evidence for the existence of a 
dominant dark matter component, but as direct evidence 
for the failure of the current Newtonian and General Rel- 
ativistic theories of gravity, in the large scale or low ac- 
celeration regimes relevant for the above situations. 

Numerous alternative theories of gravity have rece ntly 
appeared (e.g. TeVeS of a nd variations; [1], [Slf . 
F(R) theories e.g. [1^, (cj. |35|. conformal gravity the- 
ories e.g. [l9|). mostly grounded on geometrical exten- 
sions to General Relativity, and leading to laws of gravity 
which in the large scale or low acceleration regime, mimic 
the Modified Newtonian Dynamics (MOND) fitting for- 
mulas. Similarly, [221 have explored MOND not as a 
modification to Newton's second law, but as a modified 
gravitational force law in the Newtonian regime, finding 
a good agreement with observed dynamics across galactic 
scales without requiring dark matter. In fact, recently [4] 
have constructed an f{R) extension to general relativity 
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which in the low velocity limit converges to the above 
approach. 

Whilst Classical Gravity augmented by the dark mat- 
ter hypothesis provides a coherent and unified interpreta- 
tion from galactic to cosmological scales (with the inclu- 
sion of dark energy), the very profusion of modified grav- 
ity theories, mostly tested in very localised situations, 
points to the lack of any definitive theoretical contender 
to Classical Gravity. Nonetheless, a generic feature of 
all of the modified gravity schemes mentioned above is 
the appearance of an acceleration scale, ao, above which 
classical gravity is recovered, and below which the dark 
matter mimicking regime appears. This last feature re- 
sults in a general prediction; all systems where a » qq 
should appear as devoid of dark matter, and all systems 
where a << ao should appear as dark matter dominated, 
when interpreted under classical gravity. It is interesting 
that no a >> oq system has ever been detected where 
dark matter needs to be invoked, in accordance with the 
former condition. On the other hand, the latter con- 
dition furnishes a testable prediction, in relation to the 
orbits of wide binaries. For test particles in orbit around 
a IMq star, in circular orbits of radius s, the accelera- 
tion is expected to drop below qq « 1.2 x 10^^"m/s^ for 
s >7000 AU= 3.4 X 10"^pc. The above provides a test 
for the dark matter/ modified theories of gravity debate; 
the relative velocities of components of binary stars with 
large physical separations should deviate from Kepler's 
third law under the latter interpretation. 

More specifically, seen as an equivalent Newtonian 
force law, beyond s SS7000 AU the gravitational force 
should gradually switch from the classical form of = 
GM/s^ to Fmg = (GMaoy/^/s, and hence the orbital 
velocity, V'^/s = F, should no longer decrease with sepa- 
ration, but settle at a constant value, dependent only on 
the total mass of the system through V — (GA/oq)^/'*. 
That is, under modified gravity theories, binary stars 
with physical separations beyond around 7000 AU should 
exhibit "flat rotation curves" and a "TuUy-Fisher rela- 
tion" , as galactic systems in the same acceleration regime 
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do. 

An interestin g pr ecedent in this sense is given by the 
recent results of |27l| and [l^ who find evidence for a tran- 
sition in the dynamics of stars in the outer regions of a 
series of Galactic globular clusters. These authors report 
a flattening of the velocity dispersion profile in globu- 
lar clusters, outwards of the radius where accelerations 
fall below the oq threshold, in accordance with generic 
predictions of modified gravity schemes. The interpre- 
tation under Newtonian dynamics explains the observed 
flattening as due to tidal heating by the Milky Way, e.g 
[l8!], but the matter is still being debated. We also note 
the recent results of [l^ who point out various discrepan- 
cies between standard ACDM predictions and structural 
and dynamical properties of the local group, and suggest 
solutions to these in the context of modified gravity the- 
ories. 

In this paper we propose that wide binary orbits may 
be used to test Newtonian gravity in the low accelera- 
tion regime. We apply this test to the binaries of two 
very recent catalogues containing relative velocities and 
separations of wide binaries. The two catalogues are en- 
tirely independent in their approaches. The first one, [l^ 
uses data from the Hipparcos satellite to yield a moderate 
number of systems (280) relatively devoid of false posi- 
tives (10%), with a high average signal to noise ratio for 
the relative velocities of the binaries (~ 2). The second, 
[m identifies 1,250 wide binaries from the Sloan Digi- 
tal Sky Survey (SDSS) data base data release 7, which, 
compounded with a detailed galactic stellar distribution 
model, results in pairs with a very low probability of 
chance alignment (< 2%), albeit with a low average sig- 
nal to noise ratio in their relative velocities (^ 0.5). 

The paper is organised as follows: section (2) briefly 
gives the expectations for the distribution of relative ve- 
locities as a function of separation for wide binaries, un- 
der both Newtonian gravity and generically for modified 
theories of gravity. In section (3) we show the results 
of applying the test to the [2^ Hipparcos catalogue, and 
to the independent ^T] SDSS data. Our conclusions are 
summarised in section (4). 



II. EXPECTED RELATIVE VELOCITY 
DISTRIBUTIONS FOR WIDE BINARIES 

Since orbital periods for IMq binaries with separa- 
tions in the tens of AU range already extend into the 
centuries, there is no hope of testing the prediction we 
are interested in through direct orbital mapping. Fortu- 
nately, modern relative proper motion studies do reach 
binary separations upwards of 10* AU, e.g. [3, 0: [HI- 
The Newtonian prediction for the relative velocities of 
the two components of binaries having circular orbits, 
when plotted against the binary physical separation, s, 
is for a scaling of AV cx s~^/^, essentially following Ke- 
pler's third law, provided the range of masses involved 
were narrow. 



In a relative proper motion sample however, only two 
components of the relative velocity appear, as velocity 
along the line of sight to the binary leads to no proper 
motion. Thus, orbital projection plays a part, with sys- 
tems having orbital planes along the line of sight some- 
times appearing as having no relative proper motions. A 
further effect comes from any degree of orbital ellipticity 
present; it is hence clear that the trend for AV oc s~^/^ 
described above, will only provide an upper limit to the 
distribution of projected AV vs. s expected in any real 
observed sample, even if only a narrow range of masses is 
included. One should expect a range of measured values 
of projected AV at a fixed observed projected s, all ex- 
tending below the Newtonian limit, which for equal mass 
binaries in circular orbits gives: 

The problem is complicated further by the dynamical 
evolution of any population of binaries in the Galactic 
environment. Over time, the orbital parameters of wide 
binaries will evolve due to the effects of Galactic tidal 
forces. Also, dynamical encounters with other stars in 
the field will modify the range of separations and rela- 
tive velocities, specially in the case of wide binaries. To 
first order, one would expect little evolution for binaries 
tighter than the tidal limit of the problem, and the even- 
tual dissolution of wider systems. 

A very detailed study of all these points has recently 
appeared, [l6|. These authors numerically follow popula- 
tions of 50,000 IMq binaries in the Galactic environment, 
accounting for the evolution of the orbital parameters of 
each due to the cumulative effects of the Galactic tidal 
fleld at the Solar radius. Also, the effects of close and 
long range encounters with other stars in the field are 
carefully included, to yield a present day distribution of 
separations and relative velocities for an extensive popu- 
lation of wide binaries, under Newtonian Gravity. Inter- 
estingly, one of the main findings is that although little 
evolution occurs for separations below the effective tidal 
radius of the problem, calculated to be of 1.7 pc, the sit- 
uation for grater separations is much more complex than 
the simple disappearance of such pairs. 

It is found that when many wide binaries cross their 
Jacobi radius, the two components remain fairly close 
by in both coordinate and velocity space, drifting in the 
Galactic potential along very similar orbits. This means 
that in any real wide binary search a number of wide 
pairs with separations larger than their Jacobi radii will 
appear. Finally, [l^ obtain the RMS one-dimensional 
relative velocity difference, AVid, projected along an ar- 
bitrary line of sight, for the entire populations of binaries 
dynamically evolved over 10 Gyr to today, for a distribu- 
tion of initial ellipticities, as plotted against the projected 
separation on the sky for each pair. The expected Ke- 
plerian fall of AVm oc s~^/^ for separations below 1.7 
pc is obtained, followed by a slight rise in AVm as wide 
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systems cross the Jacobi radius threshold. AVid then 
settles at RMS values of « 0.1km/ s. 

These authors also explore variation in the initial (real- 
istic) distribution of semi-major axes, and the formation 
history of the binaries, finding slight differences in the re- 
sults, which however are quite robust to all the variations 
in the parameters explored, in the regime we are inter- 
ested of present day projected separations larger than 
log{s/pc) > — 2, above very small variations of less than 
0.14 in the logarithm. This represents the best currently 
available estimate of how relative velocities should scale 
with projected separations for binary stars (both bound 
and in the process of dissolving in the Galactic tides) 
under Newtonian gravity. 

The testable quantitative prediction of Classical Grav- 
ity for the distribution of data in a plot of projected AV^ 
vs. projected s is clear: one should find a spread of points 
extending below the limit defined by eq.(l}, with an RMS 
value for AVm given by the results of their figure 
(7). 

The recent proliferation of modified gravity models 
however, implies the absence of a definitive alternative to 
Classical Gravity. Further, in many cases, the complex 
formulations put forth do not lend themselves to straight- 
forward manipulations from which detailed predictions 
might be extracted for varied applications distinct from 
the particular problems under which such models are pre- 
sented. We shall therefore not attempt to test any partic- 
ular modified gravity theory, but shall only consider the 
generic predictions such theories make for "flat rotation 
curves" in the a < oq regime. That is, the predictions of 
modified gravity schemes will only be considered qualita- 
tively and generically, to first order, in terms of the upper 
envelope of observed distributions in projected AV^ vs. s 
plots to appear fiat. For circular orbits one expects: 



AVi 



MG 



C{GMao 



(2) 



where C is a model-dependent constant expected to be 
of order unity. A further correction upwards due to de- 
partures from circularity, which at this point must be 
thought of as dependent on the details of the particular 
modified gravity scheme one might pick, should also be 
included. This correction will tend to give even larger 
values of AVmg- We note that in the particular case of 
MOND, the external field effect, the fact that the over- 
all potential of the Galaxy at the solar neighbourhood 
globally puts local binaries close to the a ~ uq thresh- 
old, would imply only slight corrections om Newtonian 
predictions. 

We see that all we need is a large sample of relative 
proper motion and binary separation measurements to 
test the Newtonian prediction for the RMS values of the 
1 dimensional relative velocities of [l^, and to contrast 
the AVn oc s~^/^ and the AVmg = cte. predictions for 
the upper envelope of the AV vs. s distributions. It is 
important to have a sample as free of chance alignments 
as possible, as the inclusion of non-physical stellar pairs 



would blur the test, potentially making a conclusion sus- 
pect. Also, it is desirable to limit the range of masses of 
the stars involved, as a spread in mass will also blur any 
trends expected for the upper limit of the AV distribu- 
tions, although not terribly so, given the small powers to 
which mass appears in both predictions. 

To end this section we briefly recall the first order tidal 
limit calculation of 



dR 



Ar = 



jArf 



(3) 



which leads to the tidal density stability condition of 
Ps > p for the density of a satellite of extent Ar and 
mass Ms orbiting at a distance Rq from the centre of 
a spherical mass distribution M{R) having an average 
matter density p internal to Rq resulting in a gravita- 
tional force Fext{R), under the assumption Ar << R. 
The equivalent calculation under the force law given by 
the a « ao fimit of Fmg = (GMao)^^'^ / R is given by: 



{GM{R)aoy/^ ^ (GM,ao)'/' 



i?2 



-Ar 



leading to: 



Ps > 



Ar 



(4) 



(5) 



as the equivalent of the classical tidal density criterion, 
as a first generic approximation under modified gravity. 
Since the spatial extent of wide binaries will always be 
much smaller that their Galactocentric radii, equation (5) 
shows that under modified gravity, to first order, wide 
binaries will be much more robust to tides than under 
Newtonian gravity. In the following section we apply the 
test we have identified to two recent catalogues of wide 
binaries which became available over the previous year, 
the SLoWPoKES catalogue of SDSS wide binaries by E| 
and the Hipparcos satellite wide binaries catalogue of [29| . 



III. OBSERVED WIDE BINARY SAMPLES 
A. The Hipparcos wide binaries 

The Q catalogue of very wide binaries was con- 
structed through a full Bayesian analysis of the combined 
Hipparcos database, the new reduction of the Hipparcos 
catalogue, 32], the Tycho-2 catalogue, and the Ty- 
cho double star catalogue, [l3| mostly, amongst others. 
There, probable wide binaries are identified by assign- 
ing a probability above chance alignment to the stars 
analysed by carefully comparing to the underlying back- 
ground (and its variations) in a 5 dimensional parameter 
space of proper motions and spatial positions. The au- 
thors have taken care to account for the distortions intro- 
duced by the spherical projection on the relative proper 
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FIG. 1. The figure shows projected relative velocities and sep- 
arations for each pair of wide binaries from the [29l ] Hipparcos 
catalogue having a probability of being the result of chance 
alignment < 0.1. The average value for the signal to noise 
ratio for the sample shown is 1.7. The upper limit shows the 
flat trend expected from modifled gravity theories, at odds 
with Kepler's third law, shown by the s"^^^ solid line. 

motion measurements, A/x. When angular separations 
cease to be small, small relative physical velocities be- 
tween an associated pair of stars might result in large 
values of A/i. A correction of this effect is introduced, 
to keep Afj, values comparable across the whole binary 
separation range studied. 

We have taken this catalogue and kept only binaries 
with a probability of non-chance alignment greater than 
0.9. The wide binary search criteria used by the au- 
thors requires that the proposed binary should have no 
near neighbours; the projected separation between the 
two components is thus always many times smaller than 
the typical interstellar separation, see [l^ . We use the re- 
ported distances to the primaries, where errors are small- 
est, to calculate projected AV^ and projected s from the 
measured A^ and AO values reported by [1^. Although 
the use of Hipparcos measurements guarantees the best 
available quality in the data, we have also further pruned 
the catalogue to remove all binaries for which the final 
signal to noise ratio in the relative velocities on the plane 
of the sky was lower than 0.3. 

We plot in figure (1) a sample of 280 binaries con- 
structed as described above, having distances to the Sun 
within 6 < d < 100 in pc. The slanted line gives the 
Newtonian prediction of eq.(l) to the upper limit ex- 
pected on the relative velocities shown, which appears in 
conflict with it, as they are defined by a neat horizontal 
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FIG. 2. The figure shows projected relative velocities and 
separations for each pair of wide binaries from the 1 11] SDSS 
catalogue within the distance range (225 < d/pc < 338). The 
average value for the signal to noise ratio for the sample shown 
is 0.5. The upper limit shows the fiat trend expected from 
modified gravity theories, at odds with Kepler's third law, 
shown by the s~^^^ solid line. 



upper limit, as generically predicted by modified gravity 
theories, eq.(2). Figure (1) could then be a first direct 
evidence of the breakdown of classical gravity theories in 
the low acceleration regime oi a < uq. 

The average signal to noise ratio for the data in figure 
(2) is 1.7, with an average error on AV^ of 0.83 km/s, 
which considering a 2a factor from the top of the distri- 
bution to the real underlying upper limit for the sample, 
results in 3 km/s as our estimate of the actual physi- 
cal upper limit in AV. Comparing with eq.(2), the fac- 
tor accounting for non-circular orbits in modified gravity 
comes to 4.5. That this factor is significantly larger than 
the of Newtonian gravity is to be expected, as objects 
are much more tightly bound in MOND-type schemes. 

B. The SDSS wide binaries 

The Sloan low mass wide pairs catalogue (SLoW- 
PoKES) of 11] contains a little over 1,200 wide binaries 
with relative proper motions for each pair, distances and 
angular separations. Also, extreme care was taken to in- 
clude only physical binaries, with a full galactic popula- 
tion model used to exclude chance alignment stars using 
galactic coordinates and galactic velocities, resulting in 
an estimate of fewer than 2% of false positives. As with 
the Hipparcos sample, this last requirement yields only 
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FIG. 3. The figure shows the number of systems with AV 
values below the Newtonian upper limit prediction of eq. (1), 
as a fraction of the total per bin, for the SDSS data of 
as a function of projected binary separations. 



isolated binaries with no neighbours within many times 
the internal binary separation. We have also excluded all 
systems with white dwarfs or subdwarf primaries, where 
distance calibrations are somewhat uncertain. As was 
also done for the Hipparcos sample, all triple systems re- 
ported in the catalogue were completely excluded from 
the analysis. 

Given the large range of distances to the SDSS binaries 
(46 < d/pc < 992), we select only 1/3 of the sample 
lying within the narrow distance range (225 < d/pc < 
338), which forms the most homogeneous set in terms 
of the errors in A.V, excluding data with large errors at 
large distances. Again, we use the reported distances 
to the primaries, where errors are smallest, to calculate 
projected AV and projected s from the measured A/x, 
A6 and d values reported by to plot figure (2). The 
figure shows 417 binaries with average signal to noise 
ratio and average errors on AV of 0.5 and 11.3 km/ s, 
respectively. 

The slanted solid line gives the Newtonian prediction 
of eq.(l). It is clear that the upper envelope of the distri- 
bution of AV measurements from the catalogue does not 
comply with Kepler's third law. As was the case with 
the Hipparcos sample, the upper envelope of the distri- 
bution of observed measurements describes a flat line, as 
expected under modified gravity schemes. 

Figure (3) is a plot of the number of SLoWPoKES sys- 
tems for the full distance range, with AV values below 
the Newtonian prediction of eq.(l), and hence consistent 
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FIG. 4. The solid curve gives the RMS values for one di- 
mensional projected relative velocities as a function of pro- 
jected separations, for the detailed dynamical modelling of 
large populations of wide binaries evolving in the Galactic en- 
vironment, taken from [l^. The same quantity for the data 
from the catalogues analysed is given by the points with error 
bars; those with narrow log{s) intervals being from the Hip- 
parcos sample of [2^, and those two with wide log{s) intervals 
coming from the SDSS sample of [TH . 



with it, as a fraction of the total per bin, as a func- 
tion of projected binary separations. We see this frac- 
tion starts off being consistent with 1, but begins to de- 
crease on approaching log{s/pc) —2, after which point 
it rapidly drops, to end up consistent with on reaching 
separations of around log{s/pc) ~ —1. Again, the result 
matches the qualitative generic expectations of modified 
gravity schemes, but would call for further explanations 
under classical gravity. 

The average signal to noise values in Ay for the full 
distance range of the 11] catalogue is 0.48. The average 
error on AV for the full SDSS sample is 12 km/ s, which 
considering a 2a factor from the top of the complete dis- 
tribution to the real underlying upper limit, results in the 
same 3km /s as obtained for the [29| Hipparcos catalogue. 

We end this section with figure (4), where we calcu- 
late the RMS value of the one-dimensional relative ve- 
locity difference for both of the samples discussed, after 
binning the data into constant logarithmic intervals in 
s. This quantity is given by the points, where the error 
bars simply show the propagation of the errors on Afi 
and d, reported by the authors of the catalogues. We 
construct AVm by considering only one coordinate of 
the two available from the relative motion on the plane 
of the sky. Thus, each binary can furnish two AVm mea- 
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surements, which statisticaUy should not introduce any 
bias. Indeed, using only A/i; or only A/if, or both for 
each binary, yields the same mean values for the points 
shown. The small solid error bars result from consid- 
ering an enlarged sample where each binary contributes 
two AViD measurements, while the larger dotted ones 
come from considering each binary only once, and do not 
change if we consider only A/i; or only A/i^. The series 
of small /oq(s) interval data are for the Hipparcos cata- 
logue o f I29l . while the two broader crosses show results 
for the [Tl| SDSS sample of figure (2). For this last case, 
the much larger intrinsic errors mean that to compensate 
through weight of numbers the low signal to noise ratio 
of this catalogue, imposes the loss of separation resolu- 
tion through the use of only two bins, with vertical error 
bars which are only relevant if the sample is doubled, as 
described above. 

The solid curve is the Newtonian prediction of the full 
Galactic evolutionary model of [16] for a randomly ori- 
ented population of wide binaries with a realistic distri- 
bution of eccentricities, both bound and in the process 
of dissolving. Note that the results of this simulation 
deviate from Kepler's law for s larger than the New- 
tonian Jacobi radius of the problem of l.Tpc, whereas 
the deviation shown by the samples of binaries studied 
also occur at much smaller separations (see below) . Even 
considering the large error bars, where each binary con- 
tributes only one AVm value, we see eight points lying 
beyond la, making the probability of consistency be- 
tween this prediction and the observations of less than 
(0.272)^=3 X 10"^. The only point where this model is 
slightly at odds with the selection of the Hipparcos sam- 
ple of 29], is that j3| assume IMq stars for their bina- 
ries, while the typical mass of the stars in the binaries 
we examine is closer to 0.5Mq. This detail would only 
shift the Newtonian prediction a factor of 2^^'^ further 
away from the measurements. We obtain an RMS value 
for AViD compatible with a horizontal line at 1 km/s, 
in qualitative agreement with expectations from modi- 
fied gravity schemes. The vertical dashed line marks the 
a = ao threshold; we clearly see the data departing from 
the Newtonian prediction outwards of this line, and not 
before. 

Under Newtonian Gravity one would need to look for 
an alternative dynamical dissociation and heating mech- 
anism for the binaries we analyse which might result in 
relative velocities an order of magnitude above the re- 
sults obtained by [11] . The Hipparcos catalogue has been 
closely studied for over a decade, and not only the values 
reported, but also the uncertainties in them are now well 
established. It is highly unlikely that these confidence 
intervals might have been systematically underestimated 
by the community by a factor of 3, as would be required 
for the data in fig. (4) to be consistent with the Newtonian 
prediction of [IQ] . 

The two SDSS points are clearly consistent with the 
Hipparcos measurements, but given the much larger error 
bars, they are also marginally consistent with the New- 



tonian prediction at a 2a level. A fuller discussion of the 
robustness of our results to the various sample selection 
effects and errors involved appears in the Appendix. 

The consistency of the results obtained from both cat- 
alogues provides a check of the physical reality of the 
trends presented. The two completely independent, very 
carefully constructed catalogues, each using different sets 
of selection criteria, each perhaps subject to its own inde- 
pendent systematics, are consistent with the same result, 
a constant horizontal upper envelope for the distribution 
of relative velocities on the plane of the sky at an intrinsic 
value of 3 fcm/s±l km/s, extending over 3 orders of mag- 
nitude in s, with a constant RMS AVid value consistent 
with 1 fcm/s±0.5 km/s. This supports the interpreta- 
tion of the effect detected as the generic prediction of 
modified gravity theories. 

Given the relevance of the subject matter discussed, it 
would be highly desirable to obtain an independent con- 
firmation, ideally using a purposely designed sample in- 
cluding cuts at a variety of stellar masses, in order to test 
the scaling of the fourth power of the dynamical velocities 
with mass expected under modified gravity schemes. To 
include also radial velocity measurements would require 
sampling over many epochs, as the only way to allow 
for the effects of nearby companions in these frequently 
hierarchical systems. 

IV. CONCLUSIONS 

We have identified a critical test in the classical grav- 
ity/ modified gravity debate, using the relative velocities 
of wide binaries with separations in excess of 7000 AU, 
as these occupy the a < ao regime characteristic of mod- 
ified gravity models. We present a first application of 
this critical test using the best currently available data; 
a large sample of wide binaries from the SDSS with low 
signal to noise on the relative velocities, and a smaller 
Hipparcos satellite sample with signal to noise ~ 2 on 
the relative velocities of the binaries sampled. 

Results show constant relative RMS velocities for the 
binary stars in question, irrespective of their separation, 
in the a < ao regime sampled. This is quantitatively 
inconsistent with detailed predictions of Newtonian dy- 
namical models for large populations of binaries evolving 
in the local galactic environment. 

Our results are qualitatively in accordance with generic 
modified gravity models constructed to explain galactic 
dynamics in the absence of dark matter, where one ex- 
pects constant relative velocities for binary stars, irre- 
spective of their separation, in the a < oq regime sam- 
pled. 
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Appendix A: Calculation of confidence intervals 

We begin this section with a discussion of various anal- 
yses performed to test for the possibility that our results 
could have been driven by potential systematics and se- 
lection effects in the catalogues. 

We first test for the option that the results of figures 
(1) and (2) were distorted from the Newtonian predic- 
tion by errors which correlate with the separation of the 
binaries, increasing as the separation increases, to yield 
the trends obtained. For the Hipparcos sample we ranked 
the binaries by separation, s, and calculated the average 
errors on the resulting AV in each the tightest third, 
middle third and widest third of the binaries, yielding 
values 0.8 (0.4), 0.7 (0.11) and 1.0 (0.4) respectively, in 
km/s. The numbers in parenthesis giving the dispersion 
of the distributions of errors in each of the three thirds 
of the sample. It can be seen that there is no increase 
either in the average values of the errors in AV^, or in 
the width of the distributions of errors, with increasing 
binary separations. For the SDSS sample the average er- 
rors in AV^ show only a very slight increase of a factor of 
1.3 over the entire range probed. Thus, for both samples, 
the trends of figures (1), (2) and (4) can not be explained 
as arising from the Newtonian prediction and an increase 
in the errors with binary separation. In essence, the data 
presented are inconsistent with the Newtonian prediction 
not because of differences in the details of the trends, but 
because the former presents multiple real detections at a 
level of 3fcm/s± Ikm/s, while the latter is smaller by an 
order of magnitude. The resulting RMS values for the 
various data points are of ~ Ifcm/s ± 0.5fcm/s, while the 
Newtonian prediction lies below the detection by about 
a factor of 10. 

We next check against systematics in the catalogues, 
which would preferentially arise in the lower signal to 
noise points, or alternatively, in the most distant ones. 
For the Hipparcos sample, we repeat the experiment 
including only points with a signal to noise ratio on 
Ay > 1, and then including only the binaries nearest 
to earth with d < 50pc. These two sub- samples have 
186 and 130 binaries, with average signal to noise ratios 
on AV above the value of 1.7 of figure (1), of 2 and 2.2 
respectively, and yield results indistinguishable from figs. 
(1) and (4). Also, we checked explicitly the results of the 
SDSS sample for any systematics with distance to the 
binaries analysed, and found the flat upper envelope to 
be robust to the choice of distance range taken. 

Although no specific cut in resulting AV or even on 
measured A/x was built into the Hipparcos catalogue, the 
SDSS one includes the selection cut that the signal to 
noise ratio on A^ should be < 2. We repeat the SDSS 



analysis tracing the 15% of the points closest to this cut, 
and find that these, the ones which just made the cut, 
do not define the flat upper envelope on AV. For both 
the Hipparcos and the SDSS samples, any observational 
bias/truncation would appear in (A/i, AO) space, not in 
the {AV, s) one. 

We end with a detailed description of the calculation of 
the points in figure (4) and their error bars, directly from 
the data published in the catalogues used. Whenever 
z — fix), the error in z is given by 6z = j^jfe, where 
5x is the error in a;. If z is a function of several variables 
z — /(xi, X2, ...), the absolute error, which we take, is 
5z = I §^ which is always greater than standard 

deviation cr^ = \JYAMi^^^ ■ 

Following these rules we estimate the absolute error 
in each AV measurement, and in the RMS velocity for 
the samples studied. For the Hipparcos sample we have 
taken columns 11,12,13 and 14 of the online catalogue of 
[29| . which contains the values of the relative motions of 
each binary and their errors, as well as columns 15 and 
16, which contain the distance to the binary, d, and its 
corresponding error, 5d. 

For all calculations, we considered the distance to the 
system as the distance to the primary. For each binary 
we calculated the projected relative velocity Av in km/s 
and its error as: 



Av = 4.74 X lO^^A^d, (Al) 



5{Av) = 4.74 X lQ-^{5{Ay)d + AfiSd), (A2) 

hence, on {Av)^the error is S{Av)^ = 2AvS{Av). The 
RMS velocity is now. 
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-Y{Av)^ = v/< (At-)2 >. 



(A3) 



We then binned the data into constant logarithmic in- 
tervals in s with a width of 0.32 and calculated Avrms 
for each bin, n the number of binary systems that fall 
into each bin. The error in Avrms for each bin is now: 



S{Avrms) 



1 < SjAv)^ > 

2 ^< (Av)^ > 



(A4) 



To compare with the prediction of the RMS values 
for the one dimensional projected relative velocities of 
the evolutionary model of [l6| we construct Avjims by 
considering only one coordinate. For the case of the 
SDSS sample of [111, we have calculated the relative 
proper motion for each coordinate A/i^ = /Xq,^ — /i^^ 
and Afis = fJ'SA ~ P^Sb taking columns 18 to 25 of their 
online catalogue, which contain the proper motions and 
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their errors in equatorial coordinates for each component 
of the binary system, to estimate the projected relative 
velocity we have used column 26 containing the distance 
to the primary, in this case we consider an error of 15% in 



[1] Allen, C, Poveda, A., Herrera, M. A., 2000, A&A, 356, 
529 

[2] Angus, G.W., 2008, MNRAS, 387, 1481 

[3] Bekenstein, J., 2004, Phys. Rev. D., 70, 3509 

[4] Bernal T., Capozziello S., Hidalgo J. C, Mendoza S., 

2011, Eur. Phys. J. C, 71, 1794 
[5] Bruneton, J. P., Esposito-Farese, G., 2007, Phys. Rev. 

D., 76, 124012 

[6] Capozziello, S., Cardone, V. P., Troisi, A., 2007, MN- 
RAS, 375, 1423 
[7] Chaname, J., Gould, A., 2004, ApJ, 601, 289 
[8] Chiu, M. C., Ko, C. M., Tian, Y., Zhao, H. S., 2011, 

Phys. Rev. D., 83, 06352 
[9] Covey, K. R., et al., 2007, AJ, 134, 2399 
[10] Delfosse, X., Forveille, T., Segransan, D., Benzit, J.L., 

Udry S., Perrier, C, Mayor, M., 2000, A&A, 364, 217 
[11] Dhital, S., West, A. A., Stassun, K. G., Bochanski, J. J., 

2010, AJ, 139, 2566 
[12] Fabricius, C, H0g, E., Makarov, V. V., Mason, B. D., 

WycofT, G. L., & Urban, S. E., 2002, A&A, 384, 180 
[13] Halle, A., Zhao, H. S., Li B., 2008, ApJS, 177, 1 
[14] Hernandez, X., Mendoza, S., Suarez, T., Bcrnal, T., 

2010, A&A, 514, AlOl 
[15] H0g, E., et al. 2000, A&A, 357, 367 
[16] Jiang, Y. F., & Tremaine, S., 2010, MNRAS, 401, 977 
[17] Kroupa, P. et al., 2010, A&A, 523, 32 



the distance and we have followed the same path as in the 
previous case for calculating the RMS relative velocities 
and their errors. 



[18] Lane R. et al., 2010, MNRAS, 406, 2732 

[19] Mannheim, P. D., Kazanas, D., 1989, ApJ, 342, 635 

[20] McGaugh, S. S., Wolf, J., 2010, ApJ, 722, 248 

[21] McGaugh, S. S., 2011, Phys. Rev. Lett., 106, 121303 

[22] Mendoza, S., Hernandez, X., Hidalgo, J. C, Bernal, T., 

2011, MNRAS, 411, 226 
[23] Milgrom, M., Sanders, R. H., 2008, ApJ, 678, 131 
[24] Milgrom, M., 2010, Phys. Rev. D., 82, 043523 
[25] Sanchez-Salcedo, F. J., Hernandez, X., 2007, ApJ, 667, 

878 

[26] Sanders, R. H., McGaugh, S. S., 2002, ARA&A, 40, 263 
[27] Scarpa R., Marconi G., Gimuzzi R., Carraro G., 2007, 
A&A, 462, L9 

[28] Scarpa R., Marconi G., Carraro G., Falomo R., Villanova 

S., 2011, A&A, 525, A148 
[29] Shaya, E. J., Oiling, R. P., 2011, ApJS, 192, 2 
[30] Skordis, C, Mota, D. F., Ferreira, P. G., Boehm, C, 

2006, Phys. Rev. Lett. 96, 011301 
[31] Sobouti, Y., 2007, A&A, 464, 921 
[32] van Leeuwen, F., 2007, A&A, 474, 653 
[33] Zhao H. S., Famaey B., 2010, Phys. Rev. D, 81, 087304 
[34] Napolitano N. R., Capozziello S., Romanowsky A. 

J., Capaccioli M., Tortora C, 2012, ApJ in press, 

arXiv: 1201. 3363 
[35] Capozziello S., De Laurentis M., 2011, Phys. Rep. 509, 

167 



